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ABSTRACT 

A numerical solution to the problem of thermal elastohydrodynamic lubri- 
cation of line contacts was obtained by using a finite difference formula- 
tion. The solution procedure consists of simultaneous solution of the thermal 
Reynolds equation, the elasticity equation, and the energy equation subject to 
appropriate boundary conditions. Pressure distribution, film shape, and 
temperature distribution were obtained for fully flooded conjunctions, a 
paraffinic lubricant, and various dimensionless speed parameters while the 
dimensionless load and materials parameters were held constant. Reduction in 
the minimum film thickness due to thermal effects (as a ratio of thermal to 
isothermal minimum film thickness) is given by a simple formula as a function 
of the thermal loading parameter Q: 

^min 10 

^min,I 10 + 

Plots of pressure distribution, film shape, temperature distribution, and flow 
are shown for some representative cases. 


INTRODUCTION 

Elastohydrodynamic lubrication of nonconformal contacts has been an 
important aspect of research in the field of tribology for the last two 
decades primarily because of its paramount importance in the lubrication of 
heavily loaded contacts in machine elements, viz roller bearings, ball bear- 
ings, gears, and traction drives. The variation of viscosity with pressure 
and the elastic deformation associated with the high pressures generated in 
the contact region are the major causes of the complexity attributed to the 
problem. The first major effort in arriving at a numerical solution of the 
problem was made by Grubin and Vinogradova (1949), who incorporated both elas- 
tic deformation and viscosity-pressure characteristics in the inlet zone anal- 
ysis of the hydrodynamic lubrication of nonconformal line contacts. Elastic 
deformation in the contact region was approximated by the Hertzian deformation 
of the solid in that region. 

Dowson and Higginson (1961) obtained a numerical solution to the problem 
of elastohydrodynamic lubrication of nonconformal line contacts by incor- 
porating both elasticity and viscous effects in the coupled solution of the 
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Reynolds and elasticity equations. They obtained a formula for the minimum 
film thickness that was based on their theoretical results and that showed the 
effects of the load, speed, and materials parameters. 

Thermal effects in nonconformal line contacts were . investigated by 
Sternlicht et al. (1961) assuming that the heat generated due to viscous shear 
was carried away totally by the lubricant and that no heat was conducted into 
the contacting solids. However, in the lubrication of nonconformal contacts 
and within the practical range of loads and rolling speeds, conduction of heat 
to the solids appears to be the major mode of heat transfer rather than con- 
vection by the lubricant. Cheng (1965) presented a numerical solution to the 
thermal elastohydrodynamic lubrication of rolling-sliding line contacts that 
takes into account both the conduction of heat into the solid from the lubri- 
cant and convection by the lubricant. Results showed that temperature had no 
significant effect on the magnitude of the film thickness. However, the fric- 
tional force in the case of rolling-sliding contacts was affected by the tem- 
perature rise in the film. In a similar paper, Dowson and Whitaker (1965) 
also produced a complete thermohydrodynamic solution to the lubrication of 
rolling-sliding line contacts and concluded that the thermal effects on the 
film thickness were insignificant. However, the rolling speed and the inlet 
viscosity of the lubricant used were low, and it is not surprising that ther- 
mal effects were not observed. On the contrary, the work done by Cheng 
(1967), Murch and Wilson (1975), and Goksem and Hargreaves (1978) on the 
effect of viscous shear heating in the inlet zone in a Grubin type of analysis 
revealed a much larger influence of temperature on the film thickness, par- 
ticularly at high rolling speeds. Recently, Kaludjercic et al. (1980), in a 
rigorous simultaneous solution of the thermal Reynolds equation, the elastic- 
ity equation, the energy equation, and the transient heat conduction equation 
in the solids showed that temperature effects on the film thickness can be 
significant at higher rolling speeds and presented a formula for the reduction 
in the central and minimum film thicknesses caused by viscous heating in the 
film. 


Recently, Hamrock and Jacobson (1983) presented a solution for line con- 
tacts wherein the computations are carried out from the inlet to the outlet of 
the conjunction as one complete solution. These results formed the basis for 
the present investigation, and the numerical procedure adopted retains all of 
the essential features of the procedure developed by Hamrock and Jacobson 
(1983). The thermal aspect of the problem uses an approach similar to that 
adopted by Dowson and Whitaker (1965) or Kaludjercic et al. (1980). 

In the present work a numerical solution to the problem of thermal 
elastohydrodynamic lubrication of line contacts is obtained. It consists of 
simultaneous solution of the thermal Reynolds equation, which takes into 
account viscosity variations with pressure and temperature, the elasticity 
equation for the film shape due to a known pressure distribution, the energy 
equation for temperature due to a known pressure distribution and film shape, 
and the heat conduction equation for the solid surface temperatures due to a 
known temperature distribution in the film. Pressure distribution, film 
shape, and temperature distribution were obtained in fully flooded conjunc- 
tions for a paraffinic lubricant and for various dimensionless speed param- 
eters covering the practical range of low to moderately high rolling speeds, 
while the dimensionless load and materials parameters were held constant. In 
arriving at a converged solution for any particular case, special attention 
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was paid to assure that the mass flow rate per unit length was held constant 
throughout the conjunction. This latter effect is important, especially in 
the inlet zone of a fully flooded conjunction in which the inlet distance 
extends several Hertzian widths from the center of the contact. 

Computer plots of pressure distribution and film shape are characteristic 
of elastohydrodynamically lubricated line contacts. Thermal effects on the 
minimum film thickness are incorporated into a simple formula that shows the 
influence of the dimensionless thermal loading parameter. Sliding had no 
significant influence on the film thickness but did affect the frictional 
characteristics. 
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SYMBOLS 

2 

constant defined in eq. (26), m /N 
constant defined in eq. (37) 


b/b 


semiwidth of contact, m 

2 

constant defined in eq. (26), m /N 
specific heat of lubricant, J/kg-K 
specific heat of solid, J/kg-K 
constant defined in eq. (37) 

p 

modulus of elasticity of solid, N/m 


effective modulus of elasticity of solid, 
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dimensionless shear force per unit length 
shear force per unit length, N/m 
dimensionless materials parameter, aE' 
dimensionless film thickness, h/R^ 
dimensionless minimum film thickness 
dimensionless constant defined in eq. (13) 
film thickness, m 
minimum film thickness, m 
thermal conductivity of lubricant, W/m-K 
thermal conductivity of solid, W/m-K 
some characteristic length of contact, m 
number of nodes in semi axis of contact 
dimensionless pressure, p/E' 
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Peclet number, 
pressure, N/m^ 

2 

dimensionless thermal loading parameter, YnQU^/k 

dimensionless mass flow rate per unit length, 
mass flow rate per unit length, kg/s-m 
effective radius in direction of motion, m 
radius of curvature, m 
slide-to-roll ratio, 2(u^ - * Up) 

geometrical separation, m 
dimensionless temperature, t/tQ 
temperature, K 

inlet lubricant temperature, K 
dimensionless speed parameter, nnU./E'R 
fluid particle velocity, m/s 

dimensionless fluid particle velocity, u/ug 
solid surface velocities, m/s 
sliding velocity, u^ - u^, m/s 

surface velocity in direction of motion, (u^ + Up)/2, m/s 
fluid element volume, m^ 
dimensionless load parameter, w /E'R 
load, N/m 

dimensionless coordinate, x/b 
coordinate in direction of motion, m 

dimensionless coordinate in direction of film thickness, z/h 
coordinate in direction of film thickness, m 

2 

pressure-viscosity coefficient of lubricant, m /N 
elastic deformation, m 

temperature-viscosity coefficient of lubricant, K“^ 

1 1 

thermal expansivity of lubricant, 7 Ijt/ » K"-*- 

2 

viscosity of lubricant, N-s/m 

2 

viscosity of lubricant at ambient pressure, N-s/m 

dimensionless viscosity of lubricant, n/no 
coefficient of friction 


V Poissons ratio for solid 

3 

p density of lubricant, kg/m 

Pq density of lubricant at ambient pressure and inlet oil temperature, 

kg/m^ 

"p dimensionless density of lubricant, p/po 

3 

Pg density of solid, kg/m 

Subscripts: 
a solid a 

b solid b 

s solid 

X coordinate in direction of motion 

2 coordinate in direction of film thickness 

THEORY 

The approach used in solving the isothermal elastohydrodynamic aspect of 
the problem is similar to that used by Hamrock and Jacobson (1983). The 
approach used in solving the thermal aspect of the problem is similar to that 
adopted by Dowson and Whitaker (1965) and Kaludjercic et al. (1980). 


Reynolds Equation 

The generalized thermal Reynolds equation derived by Dowson (1962) and 
written in the present form by Fowles (1970) is written for line-contact prob- 
lems as follows: 


where 
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where 
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equation ( 1 ) can be written in dimensionless form as 
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Boundary conditions for equation (5) are 
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P = 0 ; X = 


P - 0 

^ - dr ^ 


at X = X 


exit 


(3) 

(4) 

(5) 

( 6 ) 


(7) 


Figure 1 shows the radii of the rollers used in defining equation (3). It is 
assumed that convex surfaces, as shown in figure 1 , exhibit positive curvature 
and concave surfaces, negative curvature. Therefore, if the center of curva- 
ture lies within the solid, the radius of curvature is positive; if the center 
of curvature lies outside the solid, the radius of curvature is negative. 


Film Shape 

The film shape can be written simply as 
h(x) = ho + s(x) + 6 (x) 

where 

hQ constant 

s(x) separation due to geometry of undeformed solids 
6 (x) elastic deformation 
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The separation due to the geometry of the two undeformed rollers shown in 
figure 1(a) can be described by an equivalent cylindrical solid near a plane, 
as shown in figure 1(b). The geometrical requirement is that the separation 
of the two rollers in the initial and equivalent situations should be the same 
at equal values of x. Therefore the separation due to the undeformed geo- 
metry of the two rollers can be written, using the parabolic approximation, as 

A 

Figure 2 shows a rectangular area of uniform pressure. From Timoshenko and 

Goodier (1951) the elastic deformation at a point x on the surface of a 
semi-infinite solid subjected to pressure at the point Xj^ can be written as 


_ 9 b _ 2 

6(x) = f p In (x - dx^ (10) 

If the pressure is assumed to be uniform over the rectangular area, the pres- 
sure can be put in front of the integral. Therefore equation (10) results in 
the following (Hamrock and Jacobson, 1983): 

6(x) = - PD (11) 

IT 

where 

D = b [(I - B) In (I - B)^ - (I + B) In (I + B)^ + 4B(1 - In b)] 
and 

b Hertzian semiwidth of contact 


b equals b/n 

n number of nodes within semiwidth of contact 


B equals b/b 


Now the term 6(x) in equation (11) represents the elastic deformation, at a 

point "x, caused by a rectangular area of uniform pressure and width 2b. If 
the conjunction is divided into a number of equal rectangular areas, the total 

deformation at a point x due to contributions of various rectangular areas 
of uniform pressure in the conjunction can be evaluated numerically. The 
total elastic deformation caused by the rectangular areas of uniform pressure 
within a conjunction can be written as 

«,(x)=| E J = |x-I|-1 (12) 

1 —X 9 C 9 • • « 
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equations (9) and (12) into equation (8) while writing 
the film thickness in the dimensionless form gives the following: ^ 


where 


j = 


" R " ”o ^ 


k - i + 1 


V, 


i=l,2,.. . 


(13) 


Energy Equation 

the distribution within the lubricant film Is determined from 

t energy equation subject to appropriate boundary condi- 
tions. The energy equation tor tine-contact problems can be expressed, 
follows*”^ convection across the film and conduction along the film, as 


where 
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32 


I?) - V If - If - d (t) 
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u = u. + 
0 
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[ollowing substitutions and assuming that lubricant properties such 
specific heat, and thermal expansivity do not vary 
7ess form as " "^‘^"tion can be written in dimension- 


where 


= N^„puH^4v - N, TuH^ 


aZ‘ 


cv 


3X "ac' 


iP N “ /iM V 

3X “ ^vd’^ \dZ ) 


, - PnC„R 

T't/V u = u/u^; 


^ (4 

0 \8W/ 


1/2 


ac 


e ( R ^ E ')2 


kn. 


(sw) 


1/2 


E')2 2 
— - 


vd bptQk 


(17) 


The following expressions for fluid velocity and its derivative in the dimen- 
sionless form are used in the energy equation: 
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Boundary conditions for the energy equation are given as 

T = 1.0 , X = - ; T(X,0) = Tb ; T(X,H) = Tg (20) 

where and Th are the surface temperatures of the two solids and are 
evaluated by solving transient heat conduction equations for solids for a 
known temperature distribution in the film. 


Surface Temperature Equations 


Carlsaw and Jaeger (1959) solved the transient heat conduction equation 
for the case of linear heat flow and developed the expression for the solid 
surface temperature rise of a semi-infinite body under the action of a moving 
heat source. Following this, the surface temperature rise in the thermal 
elastohydrodynamic lubrication of line contacts can be expressed as 


^b = 


(7rk_p_Ui-C_ -) 

^ s s b ps' 


1/2 


/ ■ (C 


dc 


0 (x - c) 


1/2 


for solid b (21a) 


t = 
a 




1/2 


/ ‘ 


h (x - c ) 


^ 1/2 solid a (21b) 


These equations are valid for constant thermal conductivity of the solid and 
for the dimensionless parameter 


r (V) ^ > 1“ 

In a lubricated contact and for the practical range of load and speed para- 
meters, this condition is well satisfied. 

Equations (21a) and (21b) can be expressed in dimensionless form as 
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ll/2 


X \ s ps b Si 


/■ 1 /il\ 

J H V3Z/z=i 


dc 


=0 (X - 5 ) 


1/2 


(22a) 


T^CX.H) 


a/2 


X \ s ps a Si 


f 1(^\ 
J H 


dc 


H (X -5) 


1/2 


(22b) 


Viscosity-Pressure-Temperature Relationship 

The isothermal viscosity-pressure relationship employed in the present 
analysis as proposed by Barus (1893) is given as 


n = riQ e 


,“P 


The complete viscosity relationship can be expressed as 

CaP+r(tQ-t)] 

n = tiq e 

In dimensionless form it is written as 
- . 3[GP■^,t(,(1.0-T)] 

The viscosity model used for the sliding friction calculation is 
n = exp [^ 1,0 + 0.01 GP YtQ(1.0 - T)j 


(23) 


(24) 


(25) 


Density-Pressure-Temperature Relationship 

The isothermal density-pressure relationship used by Dowson, Higginson, 
and Whitaker (1962) has been employed in this work. For the temperature vari- 
ation it is assumed that within the range of temperatures considered, a linear 
model applies. Therefore the complete density relationship used can be 
expressed as 


P = Pr 


1.0 +*^Cp^ ~ ^0^] 

In the dimensionless form it is expressed as 

p = ^1.0 + ^ CE'p) 


(26) 


(27) 
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Computational Procedure 


To solve the thermal elastohydrodynamic lubrication problem, it is neces- 
sary to solve the isothermal cases first. The computational procedure 
developed by Hamrock and Jacobson (1983) for the solution of isothermal 
elastohydrodynamic lubrication of line contacts was adopted here to arrive at 
the isothermal pressure distribution and film shape in the conjunction. The 
isothermal pressure distribution and the film shape thus obtained were used as 
initial guesses for thermal cases. Equations (5), (13), and (17) were written 
in the finite difference form by using a standard central difference scheme. 

Figure 3 shows the placement of the uniform nodes within the conjunc- 
tion. The number of nodes within a semicontact width was 120 throughout all 
calculations. However, the numbers of nodes in the inlet and outlet regions 
were different for higher rolling speed cases. The system of simultaneous 
equations ((5), (13), (17), and (22)) was solved by the standard Gauss-Seidel 
iterative method. All fluid property integrals were evaluated by using the 
trapezoidal rule. The cycle for automatic computation procedure is shown in 
the flow chart in figure 4. 

Computational difficulties were encountered in solving the energy equa- 
tion and in evaluating the surface temperature near the pressure spike. Since 
the film temperatures are mostly governed by the local viscosity and the pres- 
sure gradient, it is not surprising that such difficulties were encountered. 
Therefore, in some cases, temperatures near the pressure spike could not be 
evaluated. In general, temperatures near the spike were much higher than 
temperatures at the center of the contact. Therefore, in the computer plots 
of temperature, temperatures near the spike have been suppressed in order to 
clearly depict the temperatures before and after the pressure spike. 

In the case of highest rolling speed, the temperature calculation in the 
beginning of the inlet zone also gave rise to numerical problems. This is due 
to the very high film thickness and very high recirculation flow caused by the 
rapidly converging film thickness. Therefore, in that case, temperatures in 
the first few nodes were not evaluated. However, these numerical difficulties 
are likely to occur because of the highly nonlinear nature of the energy equa- 
tion coupled with the nonlinear thermal Reynolds equation. None of these 
difficulties affected the film thickness calculations in those cases since the 
magnitude of the film thickness and the film shape are not decided by these 
regions at all. The coupled solutions were considered to be converged when 
the combined solution conserved the flow from the inlet to the exit region 
besides satisfying the convergence criteria associated with various equa- 
tions. The numerical procedures for solving the energy equation and evalu- 
ating surface temperatures fail to converge for high rolling speeds when 
sliding is introduced, even. for slide-to-roll ratios of 0.1 and 0.2. This is 
due to high temperatures generated within the contact region and perhaps also 
to the vi scosity-temperature-pressure relationship used in the present calcu- 
lations. Dowson and Whitaker (1965) and Cheng (1965) have used a different 
pressure-temperature relationship in their solution. However, Dowson and 
Whitaker (1965) did not carry out computations for higher rolling speeds. 
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Mass Flow Rate 


The mass flow rate per unit length for the elastohydrodynamically lubri- 
cated contacts with temperature effects can be written as 


m 


" (“a “b"'3 - ™2 ^ 


(28) 


In the dimensionless form the flow rate is written as 


° »o“s''x 



(29) 


Force Components 

Normal force components per unit length acting on the solids can be 
written as (fig. 5): 

Wz = Waz = W[3z = /p dx 
In dimensionless form 


w 

W =p|- = ^/p dX (30) 

X X 

The tangential force component for solid a is zero. For solid b it can be 
expressed as 


Wbx = - /p dh 

In dimensionless form it can be expressed as 



Shear Forces 


The shear force per unit length of the roller is written as 



where 



(32) 
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Therefore the shear forces on the roller surfaces can be expressed in dimen- 
sionless form as 




dX for solid a (33) 


for solid b 


(34) 


For equilibrium to be satisfied, the following must be true: 

Fa - ^b Wfjx = 0 

*^a2 “ ''^bz = ® 

and the coefficient of friction can be written as 


u = 


fa -fb"“bx 


w 


w 


(35) 


RESULTS AND DISCUSSION 

Input data for the system of rollers, the physical material properties, 
and the lubricant physical properties used in the computer calculations are 
presented in table I. Computed isothermal film thickness, thermal film thick- 
ness, and rolling traction coefficient are shown in table II. Table III pre- 
sents the computed results of midfilm temperature rise, surface temperature 
rise, and coefficient of friction for combined rolling-sliding cases. 

Representative pressure distribution and film shape are shown in 
figure 6. Figure 6(a) shows the pressure distribution and film shape for the 
low dimensionless rolling speed of 1x10“^^, and figure 6(b) shows the cor- 
responding pressure distribution and film shape for the high dimensionless 
rolling speed of 1.344xl0~10 (or 8.0 m/sec). The pressure distribution and 
film shape exhibit the characteristic pressure spike and nip in the film shape 
near the pressure spike, where the nip is that portion of the film shape from 
the tip of the pressure spike to the outlet of the conjunction. The film 
shape in the contact region is flatter for low rolling speed than for high 
rolling speed. This is due to the greater film thickness and the pronounced 
hydrodynamic effect on the pressures in the contact region for the high 
rolling speed. It is noticeable too in the pressure distribution in the con- 
tact region, which departs significantly from the Hertzian pressures when com- 
pared with low rolling speeds. This causes a shift in the location of the 
pressure spike toward the center of the contact. However, the natures of the 
film shape and pressure profile are not altered by the thermal effects, and no 
appreciable changes are noticed when compared with respective isothermal 
solutions. 

A computer plot of the mass flow rate from the inlet to the outlet is 
shown in figure 7. Conservation of mass flow rate is evident and was used as 
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the final criterion for the converged solution. Conservation of flow in the 
inlet region is an important aspect since the pressure distribution and film 

shape in the inlet region are responsible for the film shape within the con- 

junction and to a greater extent for the minimum film thickness. 

Temperature distributions in the conjunction for pure rolling cases are 
presented in figures 8(a) and (b) for two dimensionless rolling speeds, 
lxlO"ll and 1.344xl0“10, respectively. In pure rolling cases the maximum 
temperature rise for the midfilm temperature occurs on the far left in the 
beginning of the inlet zone. Surface temperatures rise from the beginning of 
the inlet zone, reach a maximum, and then drop to a lower value before the 
start of the contact region. In pure rolling cases temperatures are much 

lower in the contact region than in the inlet zone. For higher rolling speeds 

much higher temperatures are attained, as is evident from the comparison of 
figures 8(a) and (b). Viscous shear heating in pure rolling is predominant in 
the inlet zone and is pronounced at high rolling speeds, thus resulting in 
significant thermal effects on film thickness for higher rolling speeds. 

With sliding, the temperature distribution in the conjuction is altered 
totally, as is evident from the figure 9. Sliding causes severe heating due 
to viscous effects in the contact region only, where the viscosities are high, 
and does not introduce significant heating in the inlet region. Both the mid- 
film temperature and the surface temperature rise considerably due to shear 
heating as a result of sliding in the contact region. Sliding results in 
different surface temperatures in the solids, the slower surface attaining 
higher temperature than the faster surface. However, for the sliding cases 
computed, no appreciable change was noticed in the film thickness, even for 
very high slide-to-roll ratios of 0.3 and 0.5. A decompression cooling effect 
is visible in the temperature distribution after the pressure spike. 

The effect of temperature on rolling and sliding traction is shown in 
figures 10(a) and (b), respectively. As the rolling speed and slide-to-roll 
ratio increase, traction coefficients decrease with respect to their iso- 
thermal values. 

In calculating sliding friction, however, a different pressure-viscosity 
model was used since the friction coefficients calculated with the model used 
in computing the pressure and temperature distribution were unreasonably 
high. This fact was mentioned by Dowson and Whitaker (1965) as a reason for 
using a different model for thermal analysis than for isothermal analysis. 


Minimum Film Thickness Formula 

Besides the usual dimensionless parameters associated with the isothermal 
film thickness calculations, the only other dimensionless parameters that 
influence the magnitude of film thickness in thermal elastohydrodynamic lubri- 
cation of rectangular contacts are 


(1) Thermal loading parameter 


Q = 


^^0"s 

k 


(2) Peclet number 


Pe = 


°^p“s”x 
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with the dimensionless group as Pe However, the influence of the 

dimensionless group of the Peclet number Pe on the minimum film thick- 
ness is negligible. Therefore the functional relationship can be suitably 
represented (as the ratio of thermal to isothermal minimum film thickness) by 


or 


mm 

^min,I 


1.0 

■ 1.0 + f(Q) 


^min _ 1.0 

^min,I ~ 1.0 + DQ® 


(36) 


(37) 


where the constants D and a are determined by the linear regression analy- 
sis of the computed data for the minimum thickness. 


The formula for the thermal elastohydrodynamic minimum film thickness is 
obtained as 


^min 10 

^min,I 10 + 

with a regression coefficient of 0.97811, where 

u -3 fi7 i|0*71p0.57,,-0.11 

H . T = 3.07 U G W 
mm, I 

as obtained from Hamrock and Jacobson (1983). 


(38) 


(39) 


Table IV compares the thermal and isothermal minimum film thickness 
ratios obtained from equation (38) with those calculated from the Kaludjercic 
et al. (1980) and Murch and Wilson (1975) formulas. There is good correlation 
with the values given by Kaludjercic's formula, but Murch and Wilson's formula 
overestimates the thermal effect at higher rolling speeds. It is understand- 
able since the analysis done by Murch and Wilson (1975) is only approximate. 

It is a Grubin type of analysis of the inlet zone for an incompressible fluid 
and constant surface temperature. 


CONCLUSIONS 

A numerical solution to the problem of thermal elastohydrodynamic lubri- 
cation of line contacts has been achieved. It calls for simultaneous solution 
of the thermal Reynolds equation, the elasticity equation, the energy equa- 
tion, and the transient heat conduction equation in the solids. The following 
conclusions were reached: 

1. Thermal effects on the minimum film thickness are negligible at low 
rolling speeds. 

2. At high rolling speeds there can be significant reduction in the 
minimum film thickness due to viscous shear heating in the conjunction. Ther- 
mal effects are incorporated in terms of the dimensionless thermal loading 
parameter Q into isothermal film thickness equations as follows: 
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^min 

^min,I 10 + 

3. Sliding has negligible effect on the film thickness or film shape, 

especially at low rolling speeds. . ,.j. 

4. Thermal effects considerably reduce both rolling traction and sliding 
traction, especially at high rolling speeds and for the combination of rolling 
with sliding. 
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TABLE I. - INPUT DATA, MATERIAL PROPERTIES, AND LUBRICANT PROPERTIES 


Radius of equivalent roller on plane, cm 1.11125 

Inlet temperature of lubricant, K 313 

Inlet viscosity of lubricant, Pa-s 0.0411 

Inlet density of lubricant, kg/m^ 866 

Pressure-viscosity coefficient of lubricant, a, GPa'^ 22.76 

Temperature-viscosity coefficient of lubricant, 0.04666 

Pressure-density coefficients of lubricant, GPa~^ 

A 0.582744 

C 1.68348 

Thermal expansivity of lubricant, 6.5x10”^ 

Thermal conductivity of lubricant, W/m-K 0.12 

Specific heat of lubricant, kJ/kg-K 2.0 

Thermal conductivity of steel rollers, W/m-K 52 

Specific heat of steel rollers. J/kg-K 460 

Density of steel rollers, kg/m^ 7850 

Elastic modulus of steel rollers, GPa 200 

Poissons ratio of steel rollers 0.3 


TABLE II. - ISOTHERMAL FILM THICKNESS, THERMAL FILM THICKNESS, AND 
ROLLING TRACTION COEFFICIENT 

[Dimensionless load, W, 2 . 0478 xl 0 “ 5 ; dimensionless materials parameter, G, 5000.] 


Case 

Dimensionless 

speed, 

U 

Dimensionless 
isothermal 
minimum film 
thickness, 
^min 

Dimensionless 
thermal mini- 
mum film 
thickness, 
^min 

Ratio of 
thermal to 
isothermal 
film thickness 

Isothermal 

rolling 

traction 

coefficient 

Thermal 

rolling 

traction 

coefficient 

1 

1.000x10-11 

20.015x10-6 

20.015x10-6 

1.0000 

0.537x10-3 

0.537x10-3 

2 

2,000 

35,384 

34.581 

.9773 

.683 

.683 

3 

3.000 

44.978 

43.881 

.9756 

.836 

.836 

4 

5.000 

65.144 

61.995 

.9516 

.968 

.968 

5 

7.575 

87.943 

83.123 

.9452 

1.074 

1.051 

6 

10.000 

113.660 

105.240 

.9259 

1.7097 

1.563 

7 

13.466 

138.060 

125.740 

.9108 

1.9097 

1.678 

8 

20.000 

177.880 

156.710 

,8810 

2.1610 

1.733 











TABLE III. - MIDFILH TEMPERATURE RISE. SURFACE TEMPERATURE RISE, AND 
COEFFICIENT OF FRICTION FOR ROLLING-SLIDING CASES 

[Dimensionless rolling speed, U, 1x10“^^; dimensionless materials 
parameter, G, 5000; dimensionless load parameter, W, 2.0478x10-5.] 


Slide-to- 
roll ratio 

Sliding 
coefficient 
of friction 

Dimensionless mid- 
film temperature 
rise at center of 
contact 

Dimensionless solid surface 
temperature rise at 
center of contact 

Faster surface 

Slower surface 

0 

0 

0 

0 

0 

.04 

.0145 

.528x10-3 

.124x10-3 

.124x10-3 

.08 

.0319 

3.552 

1.615 

1.638 

.16 

.0384 

13.799 

6.416 

6.736 

.20 

.0460 

20.217 

9.359 

9.975 

.30 

.0547 

38.000 

17.460 

19.320 

.50 

.0677 

1 

73.142 

33.040 

39.490 


TABLE IV. - THERMAL AND ISOTHERMAL MINIMUM FILM THICKNESS RATIOS OBTAINED FROM 
EQUATION (38) AND FROM KALUDJERCIC AND MURCH AND WILSON FORMULAS 

[Dimensionless load parameter, W, 2.0478x10-5; dimensionless materials parameter, G, 5000.] 


Dimensionless 
speed parameter, 
U 

Dimensionless 

thermal 

loading 

parameter, 

Q 

Dimensionless 

Ratio 

of 

thermal to isothermal 

film thickness 

r ec 1 e L 
group 
parameter, 
Pe w3/2 

Present 

formula 

Kaludjercic et al. 
(1980) formula 

Murch and Wilson 
(1975) formula 

1.000x10-11 

0.564x10-2 

0.883x10-2 

0.9875 


0.9940 

0.9890 

2.000 

2.256 

1.766 

.9785 


.9876 

.9764 

3.000 

5.076 

2.649 

.9705 


.9812 

.9616 

5.000 

14.100 

4.415 

.9563 


.9682 

.9299 

7.575 

32.360 

6.688 

.9401 


.9517 

.8880 

10.000 

56.400 

8.298 

.9263 


.9361 

.8489 

13.466 

102.270 

11.890 

.9083 


.9150 

.7953 

20.000 

225.600 

17.660 

.8784 


.8767 

.7040 









Hertzian zone 


Inlet 


Center of contact 


Outlet 





120 

Nodes 


120 

Nodes 


480 Nodes- 


180 Nodes 


660 Uniform nodes 


Figure 3. - Sample nodal structure for numerical calculations. 
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Figure 7. - Representative dimensionless flow and Poiseuille term for dImen 
slonless materials parameter G of 5000, dimensionless load parameter W 
of 2.0478x10 , and dimensionless speed parameter U of 1x10"^^ 
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Figure 8. - Variation of dimensionless midfilm and solid surface tem- 
perature rise along lubricant film in pure rolling for dimensionless 
materials parameter G of 5000, dimensionless load parameter W 
of 2.0478xl0‘' and two dimensionless speed parameters U. 
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Figure 9. - Variation of dimensionless midfilm and solid surface 
temperature rise along lubricant film in rolling-sliding case for 
dimensionless materials parameter G of 5000, dimensionless load 
parameter W of 2.0478x10“^ dimensionless speed parameter U of 
1x10"^^, and two slide-to-roll ratios S. 
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(a) Variation of pure rolling friction coefficient with 
dimensionless speed parameter for dimensionless 
materials parameter G of 5000 and dimensionless 
load parameter W of 2,0478xl0"^ 
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Sllde-to-roll ratio, S 

(b) Variation of sliding friction coefficient with slide- 
to-roll ratio for dimensionless materials parameter 
G of 5000, dimensionless load parameter W of 
2,0478x10 and dimensionless speed parameter 
U of IxlO'^^ 

Figure 10. - Effect of temperature of rolling and sliding 
traction. 
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